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Abstract 



By connecting realistic spin glass models at low temperature to the highly dis- 
ordered model at zero temperature, we argue that ordinary Edwards-Anderson spin 
glasses below eight dimensions have at most a single pair of physically relevant pure 
states at nonzero low temperature. Less likely scenarios that evade this conclusion 
are also discussed. 
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I. INTRODUCTION 



Rigorous [TIJEn and non-rigorous [^] analyses ruling out mean-field pictures (for an overview, 
see [[§-§[]) of short-ranged finite-dimensional spin glass models leave open the question of whether 
there is a single pair or infinitely many pure states at low temperature (assuming that spin flip 
symmetry is indeed broken). The droplet/scaling picture []7|-|9|] asserts the existence of at most a 
single pair in all dimensions [ ^ , 11 ], but the same conclusion might arise independently of the 
other predictions of droplet/scaling. On the other hand, if infinitely many pure states occur (for 
some dimension d and temperature T), they should do so in accord with the chaotic pairs picture 
P, [T2] , |r3| ], in which only a single pair is seen in typical large finite volumes, with the particular 
pair changing chaotically with volume. 

In this paper, we present an analysis that suggests that no multi-pair pictures should occur 
below eight dimensions. The analysis is based on relating non-zero temperature Gibbs states of 
Edwards-Anderson(EA) spin glass models [ ]14| ] to zero-temperature ground states of the highly 
disordered model [ |T5|JR)| ]. The number of ground state pairs in the highly disordered model is 
known to be one for d < 8 and (uncountably) infinite for d > 8 | JT3| JTV| ] . The relation is such 
that the number of ground states in the highly disordered model should serve as an upper bound 
for the number of (physically relevant) pure states of ordinary EA spin glasses at low but nonzero 
temperature. We make no claim about the number of ground states of the ordinary models at zero 
temperature (except by other arguments for d = 2 [|T8"1]). 

The idea behind the analysis is as follows. We construct a one-parameter family of EA models, 
parameterized by A(/3), where (3 is the inverse temperature, such that A — ► oo (slowly) as j3 — > oo 
in a way to be specified shortly. The construction is such that for any finite (3 > (3q, the model is 
expected to have the same thermodynamic behavior as the "ordinary" EA nearest-neighbor Ising 
spin glass in the same dimension corresponding to A(/3 ), in the sense that if the ordinary EA 
model is in a low-temperature broken-symmetry spin glass phase at (3 Q , our model should be also, 
and the number and organization of pure states in the two models should be identical. At T = 
(i.e., j3 = oo), our model becomes the highly disordered spin glass [ p~5| , |T6] ], whose ground state 
structure is known [ [T5|JT7| ]. So the behaviors of the two models may diverge at T = 0, i.e., the 
ground state structure of the highly disordered models may differ from that of the ordinary EA 
spin glass. But the crucial point is that any T > Gibbs measure, restricted to any finite volume, 
of our model becomes supported, as T — ► 0, on spin configurations that are ground states of the 
highly disordered model within that volume. 

This results in three possible conclusions, any of which is interesting. The most reasonable 
follows from the natural expectation that the number of ground states is at least as great as the 
number of pure states at low temperatures; thus the number of ground states in the highly disor- 
dered spin glass should give an upper bound to the number of low-temperature pure states in an 
ordinary EA model. (There would however be no implication for the number of ground states in 
the EA model at T = 0.) Because the highly disordered spin glass has only a single pair of ground 
states below eight dimensions [ |I5l|17yi9[ ], the conclusion would be that the EA spin glass has no 
more than a single pair of pure states for all d < 8 at very low temperature. This conclusion would 
extend to d = 8 if it were shown that the highly disordered model has only a single pair of ground 
states also in that dimension. 

We will discuss below why the analysis and this conclusion should apply only to incongruent 
[ |20| ] pure states. Two distinct pure states that are not global flips of each other may be either 
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incongruent or regionally congruent. Spin configurations chosen from incongruent pure states 
have a nonvanishing density of relative domain walls, i.e. couplings satisfied in one but not the 
other spin configuration; otherwise they are regionally congruent. If incongruence occurs, it should 
be generated using different sequences of coupling-independent boundary conditions, but to see 
regional congruence should require choices of boundary conditions that are coupling-dependent. 

Two other possibilities, less straightforward but also intriguing, could conceivably occur. One 
of course is that the supposition that the number of (incongruent) ground states at T = is no 
smaller than the number of (incongruent) pure states at low T > does not hold for this class of 
models. Another is that there may be no continuity or even monotonicity at all in the number of 
pure states as a function of temperature in the spin glass — perhaps the most interesting but also 
least likely of the three possibilities raised here [^T|]. 

All three of these conclusions will be discussed in greater detail; we now turn to a discussion 
of the construction of the class of models to be analyzed. 



H. THE MODEL 

We will study the Edwards- Anderson (EA) Hamiltonian [[14]] on Z d , the cubic lattice in d 



dimensions: 

= - JxyCxVy , (1) 

where J denotes a realization of the couplings J xy and where the brackets indicate that the sum 
is over nearest-neighbor pairs only, with the sites x, y E Z d . We will take the spins a x to be Ising, 
i.e., o x = ±1. 

Eq. (|l]) is the EA Ising Hamiltonian for an infinite-volume spin glass on Z d ; we also need to 
define the EA model on a finite volume, given specified boundary conditions. Let be a cube 
of side 2L + 1 centered at the origin; i.e., A L = {—L, —L + 1, • • • , L} d . The finite-volume EA 
Hamiltonian is then just that of Eq. ([I]) confined to the volume A L , with the spins on the boundary 
8Al of the cube obeying the specified boundary condition. (The boundary 8Al of the volume A^ 
consists of all sites not in A^ with one nearest neighbor belonging to A^.) 

The couplings J xy are quenched, independent, identically distributed random variables; 
throughout the paper we will assume their common distribution to be symmetric about zero. Most 
studies use either the Gaussian or ± J distributions, under the assumption that the qualitative ther- 
modynamic properties in fixed dimension — existence of a phase transition at some T c (d) (whose 
value, but presumably not its existence, will depend on the nature of the distribution), the pre- 
sumed broken spin-flip symmetry (i.e., a nonzero EA order parameter qEA) below T c , the number 
of pure states below T c , and so on — will be the same for any "reasonable" coupling distribution. 
Because this remains an assumption, there is no precise definition of "reasonable", but the expecta- 
tion is that any distribution that is symmetric about zero and falls off sufficiently quickly for larger 
coupling magnitudes will all exhibit similar spin glass behavior. Thus, for example, a uniform 
distribution supported on [—1, 1] is expected to have the same basic thermodynamic properties as 
the ± J or Gaussian spin glasses. We proceed using this assumption (or weakened versions of it). 

We depart from previous studies in parametrizing the coupling distribution through the variable 
X(/3), with A — ► oo as (3 — ► oo. The couplings J xy are defined through the relation (cf. [|T5|]) 
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JT = e xy c m (2) 

where cx(/3) is chosen to ensure that the model has a sensible thermodynamic limit (i.e., a finite 
energy/spin) when [3 — > oo, and e xy and K xy are two sets of independent, identically distributed 
(i.i.d.) random variables. Each variable e xy takes on one of the two values ±1 with equal proba- 
bility, and the K xy can be chosen from any continuous distribution (e.g. uniformly from [0, 1]) so 
that the distribution of the J xy s is "reasonable", as discussed above. Two examples will be given 
shortly. 

Our general approach is to start with a "typical" EA model at A(/3o) = 1, and then embed this 
within a one-parameter family of models that has a known ground state structure (corresponding to 
A(oo) = oo). Let us denote by J a given random coupling at the origin connecting to a specified 
nearest-neighbor site. By Eq. (@), for a given realization of the random variables e and K, we have 
I = e ci e~ K ° (ci can be taken equal to 1). After a rescaling of the temperature, we 

have 

J A = e c A |J A |= eo c A e- AA ' . (3) 



We will see that a choice of c\ ~ 1/| Jq\ x , where an overbar denotes an average over coupling 
realizations, will ensure a sensible thermodynamics for the model in the limit of zero temperature 
(providing A increases slowly as (3 — > oo). 

Note that when A = the model becomes the ± J model for any starting choice of coupling 
distribution. Of more relevance is the opposite limit: as A — > oo, all initial coupling distributions 
merge into the highly disordered model discussed in |JT3| — |T7|] . In this limit the coupling distribution 
is infinitely "stretched" nonlinearly, so that every coupling magnitude occurs on its own scale. A 
full analysis of the ground state structure of this model, and the transition in ground state pair 



multiplicity, is given in QT5|JT7| ]. A quantum version of the highly disordered limit has been used 
to study the random quantum Ising model in a transverse field P2"[]. 

So given a wide range of choices for the distribution of the K's, the spin glass model described 
here has a "reasonable " coupling distribution for any (3 < oo; we will in fact work with coupling 
distributions that have a large-magnitude cutoff at finite values (depending on (3), which guarantees 
this. As (3 increases, the distribution becomes increasingly stretched, but retains its finite large- 
magnitude cutoff for any (3 < oo, so that any of these models retains the qualitative thermodynamic 
behavior of the ordinary EA model. 

As we will see, the dependence of A on (3 can be chosen (with A increasing sufficiently slowly) 
so that if an ordinary EA model, corresponding to some finite A (e.g., A(/3 ) = 1), is in a low- 
temperature spin glass phase, then so is the model described above with running A. We will 
further see that if A increases slowly enough with [3 (depending on the choice of distribution for 
K), then the Gibbs measures in the limit (3 — > oo are supported only on spin configurations that 
are ground states of the highly disordered model. That is natural since taking A(/5) — > oo slowly 
enough is roughly equivalent to first taking the (3 — > oo (T — ► 0) limit and then the A — > oo limit. 



A. Two examples 

In this section we give two examples of coupling realizations that can be used in subsequent 
analyses. As already noted, the precise form of the distribution is unimportant as long as the 
requirements listed earlier are met; all such models should exhibit the same positive-temperature 
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behavior and will merge at zero temperature. The first example starts with a flat distribution for 
J xy in the interval [—1, 1] at A(/3 ) = 1. We need to determine the distribution for K necessary to 
recover this flat distribution for J. Let 



U x = e- XK (4) 

with U 1 = e~ K having a flat distribution in [0, 1]. 

It follows that K is taken from an exponential distribution on (0, oo), with P(K) dK = 
e~ K dK. Now let J xy = e xy c x U xy ; then M — \J xy \ has the distribution 

P(M) dM = I ^(M/c^M- 1 dM, if < M < c x , (5) 
1 0, otherwise . 



It remains to choose c\. The procedure for doing this will be discussed in Sec. |IV B| ; we here give 

tA _ \ 7 1\ 
'xy ~ A ^xy L^xy 



the result, which is c\ = X, so that J x = X e xy U x and 



P(M)dM={( 1 A)[M/A] 1/A M" 1 dM, ifO<M<A, (6) 
[ 0, otherwise . 

The second example starts with a uniform distribution for K on [0, 1]. Again 

^xy A £ X y 6 y , (7) 

where we have already chosen c A = A, but the distribution of M = \J xy \ is now given by 

P(M) dM = I (VA)M- 1 dM, if Ae^ < M < X , 
I 0, otherwise . 

The coupling distribution for this case is graphed in Fig. 1. 
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FIG. 1. Sketch of the coupling density P(J xy ) dJ xy from Eq. (g) for fixed A < oo, corresponding to a 
uniform distribution of K xy on [0,1]. 

It is important to note that as A — > oo, the distributions of J* in the two examples approach 
each other; this is true in general. Both of these distributions have a temperature-dependent cutoff 
at large coupling magnitudes. While distributions with cutoffs are chosen for convenience, one 
could also choose a Gaussian distribution for J xy at A(/3 ) = 1. Because all of these distribu- 
tions are expected to give the same thermodynamic behavior, however, we will henceforth use the 
simplest of these, given by Eqs. (0) and (||). 

IE. PURE STATE STRUCTURE WITHIN THE METASTATE 

We will study this class of models and will show that (for slow enough increase of A with (J), as 
(3 — > oo, the Gibbs measures (restricted to arbitrary fixed volumes A^) become supported on the 
ground states of the highly disordered spin glass. We require that c\ be chosen so that the model 
will be sensible in the zero-temperature limit, in that that the energy density converges to a finite 
nonzero value, and further that the choice of A(/3) be made in such a way that, for each (3 > (3q, 
the model always remains within the low temperature broken-spin-flip-symmetric spin glass phase 
(assuming the corresponding ordinary EA model at (3q is in such a phase). 

If these requirements can be met, then the most natural conclusion is that in any dimension d, 
the number of ground states of the highly disordered model provides an upper bound to the number 
of pure states observed at low (but nonzero) temperature in the ordinary EA spin glass. (We 
emphasize again, though, that other conclusions remain possible, though perhaps less plausible; 
these will be discussed in Sec. |VT|.) 

As we have emphasized in earlier papers, if there are multiple pure states, the interesting, 
and physically relevant, situation is the occurrence of incongruent states. Regional congruence 
is of mathematical interest, but to see it would require a choice of boundary conditions carefully 
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conditioned on the coupling realization J . It is not currently known how to choose such b.c.'s 
in spin glasses. Numerical treatments that look for multiple pure states implicitly search for in- 
congruent ones. (It is interesting to note that recent numerical studies [^-^6]], some employing 
coupling-dependent bulk terms in the Hamiltonian, have suggested the possibility of observation 
of regionally congruent states; but see also [ [27j ] for a different interpretation of the numerical data. 
However, recent work indicates that the energetics of the interfaces found in these studies are 
inconsistent with regionally congruent pure or ground states 

We expect our analysis to hold for the number of (incongruent) pure states within any of the 
coupling-independent boundary condition metastates. The concept of metastate was introduced 
and discussed in Refs. [[Zp and and shown to be equivalent to an earlier (but somewhat differ- 
ent) construct in Ref. [[23l]. 

A metastate is a measure on (infinite-volume) Gibbs states at fixed temperature that is con- 
structed via an infinite sequence of volumes A L , with specified boundary conditions on each dA L 
chosen in a coupling-independent manner. Roughly speaking, the metastate provides the probabil- 
ity (for varying large L) of various Gibbs states appearing within A L . The equivalence of metas- 
tates constructed with certain different coupling-independent boundary conditions was shown in 
Ref. ®. 

If there are infinitely many (incongruent) pure states, a metastate should be dispersed over 
them, giving their relative likelihood of appearance in typical large volumes. If there is no incon- 
gruence, the metastate should be unique and supported on a single pure state pair, and that pair 
will appear in most (i.e., a fraction one) of the A^'s. Regionally congruent states, if present, would 
be "invisible" in the metastate, i.e., would appear in a vanishing fraction of all the A^'s. Hence, 
even if regional congruence existed, there would still be two special or "preferred" pure states that 
would be seen in a typical A L with L large. 

An example of regional congruence is afforded by interface states in ferromagnets, which can 
only be seen with carefully chosen boundary conditions, such as those of Dobrushin [j29[]. These 
are b.c.'s in which the boundary spins above the "equator" (a plane or hyperplane parallel to two 
opposing faces of A^ and cutting it essentially in half by passing at a finite height above or below 
the origin) are chosen to be plus and the boundary spins below the equator are minus. Here the 
special pair of states consists of the uniformly magnetized up and down states, which are seen in 
large volumes with random boundary conditions as well, of course, as with periodic or free b.c.'s. 
(They would also be seen even with antiperiodic b.c.'s, since the interface would pass through any 
fixed finite region only in a vanishing fraction of volumes.) An important difference with the spin 
glass case is that in the latter there is no known procedure for obtaining boundary conditions that 
would "see" any analogous regionally congruent states. 

In the current context, one needs to specify which Gibbs measures are being examined at 
T > 0. The procedure in which one chooses coupling-independent boundary conditions and then 
varies the temperature is well-defined when one is studying properties of the metastate as a whole. 
As noted, this will provide information on the number of incongruent pure states (in the metastate) 
at a given temperature. If one wants to push further in order to study the possible existence of 
regionally congruent pure states, then one needs to pick out these Gibbs states through a choice 
of coupling-dependent boundary conditions. In order to take the (3 — > oo limit, one would then 
have to change these coupling-dependent boundary conditions in some unknown way; in this case, 
the procedure of taking the zero-temperature limit is a priori not well-defined in general. For spin 
glasses, the only well-defined procedures known at this time use coupling-independent b.c.'s, as 
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T -> 0. 

The consequence is that this argument is best able to provide information on the number of in- 
congruent states in the spin glass metastate at various positive temperatures, but that it is harder to 
draw conclusions on the possible existence or nonexistence of regionally congruent states. These 
latter are the interface states that are "invisible" in coupling-independent metastates, and in any 
case, as we have argued in earlier papers [|2|j3|jr8|], are unlikely to be physically relevant even if they 
exist, since boundary conditions used in (either numerical simulations or) laboratory experiments 
on spin glasses are coupling-independent — i.e., they are not tailored to the microscopic disor- 
der configuration. (For other results and discussion concerning the distinction between coupling- 
dependent and independent boundary conditions, see [BOB]]].) 



IV. ANALYSIS 

We now analyze the behavior of our class of models both at large finite A (i.e., at low temper- 
ature) and in the A — > oo limit. As noted above, for specificity we will study the behavior of the 
model with coupling distribution given by Eqs. ([7]) and 

Our argument relies on answering three questions (in the affirmative). For a given model 
(i.e., specified coupling distribution) can one choose a c x and a A(/3) (independently of particular 
coupling realization) so that: 

1) the model has a thermodynamically sensible scaling limit as (3 — ► oo, in that the energy per spin 
converges to some C as (3 — > oo, with < C < oo; 

2) A(/3) scales in such a way to ensure that, if the model is in its low temperature thermodynamic 
phase originally, i.e., at A(/3 ) = 1, the model at rescaled (lower) temperature remains in that same 
phase; 

3) A(/3) increases slowly enough with (3 to ensure that any Gibbs state at positive temperature is 
increasingly supported, as temperature is lowered, on spin configurations that are ground states of 
the highly disordered model? 

In the following three subsections, we provide an analysis that answers these questions in the 
affirmative, and shows how such a c x and A(/3) can be chosen. 

A. Effective temperature and couplings 

In this subsection we present a simple transformation on the (inverse) temperature and the 
couplings that enables us to map the class of models under study onto more familiar ones. In the 
absence of a scaling factor c\, the Gibbs weighting factor at inverse temperature f3' is 



exp 



(3 ^ ^xy & V @x < 3y 



-\K X 

(xy) 



(9) 



We transform the temperature and couplings using 

& e xy e- A <0*- = P JJW (10) 

where (3' = /3c\(p) and 



8 



J^ = e xy c m e-*W>*«, 



(ID 



so at a given the model maps onto an equivalent one at effective temperature (3 and with effective 
couplings J£ff>. Although simple, this rewriting of the Gibbs factor provides a natural separation 
of the various factors that allows us to choose the prefactor c\ in a coupling-independent way, so 
that the energy density has a sensible f3 — > oo limit. 

To achieve the correct scaling limit, we need to choose c\(p\ so that it scales (with (3 or (3') as 
the inverse of (minus) the energy density e(/3', A) of the model with Gibbs factor (Q). This energy 
density e( / 9 / , A) is the same for all (infinite- volume) Gibbs states and for almost all realizations of 
the couplings (under the assumption that \J£ y \ < oo) and so equals its disorder average e(/3', A). 
Since (3' (and (3) will be chosen to scale to infinity rapidly as functions of A, we will choose c\ 
to scale like the inverse of (minus) e(oo, A), the ground state energy density, which in turn scales 
like (minus) | J£ y \ as A — > oo. To justify this choice of C\, we need to compute e(/3', A) in the 
low-temperature limit. This is done in the next section. 



B. Thermodynamic behavior of the zero temperature limit 



In this subsection we estimate the disorder-averaged energy per spin e((3', A). Although (3' and 
A will later be taken to depend on each other, for the purposes of this section we treat them as 
independent variables. One of the purposes of this calculation is to provide information on how 
slowly A needs to vary with j3' (and hence with (3) in order to have sensible thermodynamics in the 
zero temperature limit. In Sec. |IV D| we will study how A needs to vary with (3 in order for pure 
states at positive temperature to be supported on ground states of the highly disordered model. We 
will see that if A grows slowly with (3' and (3, these two calculations result in mutually compatible 
ranges for allowed scaling behaviors. 

We now proceed to show that the energy per spin, at large (3, can be computed as a disorder 
average over a single (arbitrary) coupling Jq = e~ XK °, which as before denotes the coupling 
(in J') connecting the spin at the origin with one of its nearest neighbors. Let Pp^x denote the 
probability (within a Gibbs measure at fixed (3'). Then, in dimension d, 

- e(/3', A) = d (Pp>,\(J$ is satisfied) - P p >, x {Jq is unsatisfied)) | J A | . (12) 

At large A (and larger (3'), the main contribution to the energy density arises from couplings 
corresponding to the smallest K X y§, and their probability of being satisfied approaches one. More- 
over, 

= - e (0, A) < -e(/3', A) < -e(oo,A) < d~\Jf\ , (13) 

where the last inequality is because not every coupling is satisfied at zero temperature. So it will be 
sufficient for our purposes to derive a lower bound on — e(/3', A) that approaches d \ Jq | as (3 — > oo. 
In arriving at a lower bound, the following general inequality will be useful: 

Pp',\(Jo is satisfied) — P^^\{Jq is unsatisfied) 

= 1 — 2Pp> ! x(J x is unsatisfied) 

> i n p 0'A J o is unsatisfied) 

— Pp'.xiJo is satisfied) ' ^ ' 
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We use the coupling magnitude distribution Eq. (Q), which arises from the flat distribution for 
K on the unit interval [0, 1]. Because most of the contribution to the energy density comes from 
couplings corresponding to small K X y's, we break up the calculation into three cases according to 
the realization of K and its 2(2d — 1) adjacent couplings: 
Case I: K > 5, where 5 is fixed and < 8 <C 1; 

Case II: K < 5, and one or more adjacent couplings have magnitudes corresponding to K < boS, 
where 6 is a constant with 1/5 > bo > 1; 

Case III: K < 5, and all adjacent couplings have magnitudes corresponding to K > b 5. 

If we denote by ej the contribution to (minus \jd times) the energy density from coupling 
realizations corresponding to Case I, and similarly for the other cases, then trivially 



- (l/d)e(/3', A) = ej + e n + e IH . (15) 

We now study the three cases separately, denoting by X/, X H and X in the indicator functions 
that equal 1 (otherwise 0) only for those coupling realizations satisfying respectively the require- 
ments of the three different cases. 
Case I: Here we have 



e/<X 7 |J A | = / e- XK dK = 0((l/X)e- xs ) (16) 



as A — > oo. 

Case II: The calculation here is similar; we have 



auuii ncie is simiiai, we nave 

l H \Jo\ < 2(2d - l)b 8 [ S e~ XK dK = 0(5) f e 

Jo Jo 



en < Iii\Jo\ < 2(2rf - l)b 5 / e~ XK dK = 0(5) / e~ XK dK (17) 

In In 

as A — ► oo. 

Case III: This case is more involved. There are four possible configurations for the two spins 

coupled through Jq ; two of these correspond to Jq satisfied and two to Jq unsatisfied. Consider 

t • Pri \Un is unsatisfied) . . _ .n-^ ™ ■ ■ ■ ■ • , i i ? ti 

the ratio p - x ■ ■ s < that appears in Eq. (|14[). This ratio is maximized by the following 

"worst case" scenario: in changing a satisfied configuration (for Jq) to an unsatisfied one, all of 
the 2d — 1 adjacent couplings touching the flipped spin change from unsatisfied to satisfied. This 
case maximizes the cost in the Gibbs factor for the change in spin configurations. It follows that 
under the requirements of case III, 



P/3',x(Jq is unsatisfied) 
Pp',x{Jo is satisfied) 



< expl-2p>e~ x5+ ( 2d - 1 W e - b ° XS ] . (18) 



Therefore 



e HI = X in | Jq |(1 — 2 P/3' : \(Jq is unsatisfied) ) 

= [1 - 0(5)] (l - Oie-W'^+W-W^ ™}) J 5 e- XK dK (19) 

as A — ► oo. 

The behavior of the factor e _2/3 e +(2d-i)2/3 e o ^ at a pp ears m gq ( Jf^ as /5' ^ oo and 
A — >• oo is sensitive to the dependence of these two parameters on each other. We have 
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e 



-2f3'e' x5 +(2d-l)2/3'e- b xs _ e -2/3' e " A ' 5 (l-(2d-l) e -( !, o-i)A 1 5 ) ^ e -2/3'e~ xs as \ ^ ^ (20) 



regardless of the detailed behavior of X((3). Furthermore, [3'e~ x& — ► oo for any 5, if A — > oo 
slower than log(/3') (i.e., if A/log(/3') — > 0). We will therefore require (for this particular coupling 
distribution, but the calculation is similar for others) that 

X = o[\og((3)} as/3^oo, (21) 

which, as we shall see, will also guarantee that A = o [log(/3')]. 

Returning to the comparison of the energy densities (as (3' — > oo) for the three cases, we see 
that ej is reduced from em by a factor of order e~ xs , and ejj is reduced from em by a factor of 
order 5. We therefore find that 

- eW, A) = d[l ± 0(5)] e~ AX dK , (22) 

so that in the joint limit (3' — > oo, A — > oo (and with the condition A = o [log(/3')] satisfied for if 
uniformly distributed on [0, 1]), since 5 can be chosen arbitrarily small, it follows that 



- e(P',\) / (d\4\) -> 1 as/3' -> oo . (23) 
Therefore, if K G [0,1] uniformly, we have 

- e((3', A) ~ |JnN = e 11 ^ = C e~ XK dK ~ 1/A as A -> oo , (24) 

Jo 

and we set c\ = A for this distribution, as in Eq. (0). Thus (3' = (3c\ = (3\ and (^TJ) will 
indeed guarantee A = o [log (/?')]. A similar calculation for the distribution of Eq. discussed 
in Sec. P A|, results in c\ ~ A + 1, and so we may choose c\ = A. More generally, one chooses 
c\ ~ 1/| Jo |, as discussed in Sec. 0, to obtain a finite, nonzero energy density as temperature goes 
to zero. 



C. Comparison to phase behavior of ordinary EA models 



We have shown that the first of our requirements, that the thermodynamics of our class of mod- 
els behave properly in the zero temperature limit, can be met, and have shown how to compute 
the energy density at low temperatures. We turn now to our second requirement, namely that the 
model remain in the low temperature spin glass phase as T is lowered from a starting temperature 
already within the spin glass phase for the X((3 ) model. Of course, the low-temperature behavior 
of the ordinary EA spin glass in dimensions three and higher is not well understood, and it could 
conceivably be the case that in some or all dimensions between three and eight it undergoes a 
succession of phase transitions (either at discrete temperatures or continuously), or has no phase 
transition at positive temperature at all (i.e., remains paramagnetic down to zero temperature), or 
has other, perhaps more exotic, behavior. Our only goal here is to show that our class of models 
behaves similarly to the ordinary EA spin glass at least for very low, nonzero temperatures. To 
keep matters simple, we'll assume that for some range of d, the ordinary EA spin glass under- 
goes a phase transition at some T c (d) > 0, such that below that temperature there is a spin glass 
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phase with broken spin flip symmetry, and with the number (and organization) of pure states not 
depending on T, for < T < T c {d). 

That our class of models should behave similarly (at (3 < oo) to ordinary EA models, in terms 
of numbers and organization of pure states as (3 changes, is not immediately obvious because the 
couplings in our models are temperature-dependent. It does seem reasonable to expect though 
that, so long as the couplings depend weakly on f3, our models should behave similarly to more 
conventional ones. However, we can improve on this expectation and show that in fact this follows 
from a natural universality hypothesis. 

The universality we have in mind is that the above assumption about ordinary EA spin glasses 
is valid separately for each fixed finite A. Thus there will be a critical inverse temperature /3 C (A) < 
oo for the Hamiltonian Eq. (P with couplings J xy = e xy c\e~ XKxy . In order to satisfy our second 
requirement, A has to grow slowly enough with (3 (or equivalently, (3 has to grow rapidly enough 
with A) so that for (3 > fi , 

/3(A) > (3 C (X) . (25) 

So, for example, if A(/3 ) = 1 and the corresponding spin glass model is in its low T spin glass 
phase, the inequality ( |2"5j ) will guarantee that the sequence of models with running A > 1 remain 
within their corresponding low T spin glass phase. The above inequality can always be satisfied 
consistently with the other conditions 1) and 3); all that needs to be done is to choose A to grow 
sufficiently slowly with (3. The remaining question will be whether the constraints imposed by 
conditions 1) and 3) are compatible; we reserve that for the following section. 



D. Ground state behavior of zero temperature limit 

We now turn to a central question in the analysis of our models, which is how slowly A must 
vary with (3 in order that condition 3) above is satisfied, i.e., whether the Gibbs states of our 
models at positive temperature are increasingly supported, as temperature is lowered, on spin 
configurations that are ground states of the highly disordered model, and whether this constraint 
is compatible with condition 1). 

Consider again a fixed (arbitrary) volume centered at the origin, with any fixed boundary 
condition chosen independently of the coupling realization. Let AE(X) denote the energy differ- 
ence between the lowest-energy state in A L at given A and the first excited state. Then condition 
3) is satisfied if (3AE{X{(3)) — > oo as (3 — > oo, because the lowest-energy state goes to the ground 
state of the highly disordered model (for that b.c, in that volume) as A — > oo. 

We again study this question in the context of a particular model, namely where K is chosen 
uniformly from the interval [0, 1]. We have here that J x ~ A e xy e~ XKxy . As A gets large for fixed 
L, the first excited state corresponds to the (multi-)spin flip that changes only the smallest magni- 
tude coupling in the (wired b.c.) invasion tree (see Refs. JT5|JT7| ]) from satisfied to unsatisfied. But 
for the distribution chosen, the magnitude of this coupling is larger than A exp[— A], so that at any 
A, AE(X) > Xe- X . Therefore, 

(3AE(X) > (3Xe~ x , (26) 
so if A grows as log(/3) or slower, condition 3) is satisfied (for this particular distribution of K). 
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Condition 1) requires a slightly stronger constraint, given by ([21]), which is certainly compat- 
ible with the constraint given above. It is gratifying that both conditions require A to grow slowly 
with (3, as initially anticipated. It might be that condition 2) would require A to grow even more 
slowly with (3, but in any event poses no conflicts with the other constraints. 

The procedure given here can be adapted for other distributions, but it is already sufficient that 
an explicit example can be constructed of a model with all of the desired properties listed earlier. 



V. DISORDERED FERROMAGNETS 

The basic argument of this paper is that for the models with temperature dependent couplings 
as in Eq. @ (and with properly chosen c\ and A(/3)), the absence of more than a single pair of 
ground states at (3 = oo for d < 8 is evidence for no more than two pure states in ordinary EA 
spin glass models with fixed couplings at very low temperatures. 

A potential flaw in this line of reasoning seems to arise in the case of disordered ferromagnets, 
corresponding to the elimination of the random signs e xy in Eq. @. After all, highly disordered 
ferromagnets also have only a single pair of ground states for d < 8 [|T5],|T7]], but ordinary dis- 
ordered ferromagnets are expected to have multiple pure states for d > 4 []32p — namely, the 
interface states [ |3"3| , |3"3] ] obtained by using Dobrushin boundary conditions. Further, this expecta- 
tion is strongly supported by the rigorous results of Bovier and Kiilske concerning the existence 



for d > 4 of such interface states in SOS models (p5Q. Indeed, the existence of such interface 
states in ferromagnets motivated Bovier and Frohlich to argue that there are more spin glass pure 
states in d > 4 than d = 3 (see Sec. 6.4 of [|3~6l]). So we need to ask why our arguments, when 
applied to disordered ferromagnets, do not lead to a contradiction for 4 < d < 8. 

In fact, they do not, for at least two reasons. One of them, already discussed above, is the lack 
of appearance of interface states in the metastate (with coupling independent boundary conditions). 
We have not argued, even for spin glasses, that there are absolutely no more than two pure states 
for d < 8, but only that not more than two appear in the metastate (and that these are the physically 
relevant ones). For both ordinary EA spin glasses in d < 8 and ordinary disordered ferromagnets 
in all d, the metastate with free or periodic (or random) b.c.'s should exhibit no more than a single 
pair of pure states. Indeed, in the highly disordered ferromagnet even for d > 8, the free or 
periodic (but not random) b.c. metastate is supported only on the single pair of homogeneously 
magnetized ground states [|T3],[F7]] . 



But one might object that in the ferromagnet setting, unlike the spin glass case, one does know 
how to choose b.c.'s, namely of the Dobrushin type, so that one can obtain Gibbs states with 
interfaces in the thermodynamic limit and metastates supported on them. For example, by taking 
a mixed b.c. with the Dobrushin "equator" between the plus and minus parts of the boundary 
occuring at a variety of heights with various weights, one could obtain a metastate supported on 
many different pure interface states. So, wouldn't the fact that for 4 < d < 8, the zero temperature 
limit of the metastate could not be supported on more than two ground states, still contradict our 
reasoning that the number of ground states in the metastate at zero temperature is an upper bound 
for the number of pure states at very low temperature? 

Here is where the second reason for a distinction between ferromagnets and spin glasses comes 
into play. It concerns the universality hypothesis of our argument, that for any fixed A < oo, the 
low temperature phase structure of a spin glass should be qualitatively the same, regardless of 
the value of A. For spin glasses, this seems a perfectly plausible working hypothesis; indeed, if 
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this were not so, then, e.g., a Gaussian distribution could have a different low T thermodynamic 
structure, in terms of pure state multiplicity, than a uniform distribution on [— J, J], or a ±J spin 
glass. But this is not likely; it is expected that only non-universal features such as T c (d) depend on 
the coupling distribution (assuming that it is symmetric about zero and without slowly-decaying 
tails). 

But this is not so for disordered ferromagnets, where there are no energy cancellations along an 
interface — unlike in spin glasses. This may well lead to a greater sensitivity of interface stability 
to the strength of the disorder. It seems quite plausible, as suggested to us by Bovier and Kiilske, 
that for some dimensions d > 4, the usual interface states obtained through normal Dobrushin 
b.c.'s could perhaps disappear for A above some critical value A c < oo, even for arbitrarily low 
temperature. That is, for A < A c , the interface would be flat for both T = and small T > 0, but 
for A > A c , the interface would be rough for any T > (and perhaps also for T = 0). 

Our analysis suggests that this is indeed so for d < 8, so that there would be only two pure 
states (even for the mixed Dobrushin b.c. metastate described above) for ordinary disordered ferro- 
magnets at very low temperatures, providing the disorder is sufficiently (but not infinitely) strong. 
For 4 < d < 8, this could be a consequence of the strong disorder either completely destabilizing 
the interface states so that they are entirely absent, or else of partially destabilizing them so that 
finding them would require Dobrushin type b.c.'s, but with a coupling-dependent equator (at a 
nonconstant height). 

VI. DISCUSSION 

In this paper we have constructed a class of nearest-neighbor spin glass models with (unusual) 
coupling distributions depending on a disorder strength parameter A that itself is temperature de- 
pendent. These models are designed to have the property that their thermodynamic equilibrium 
properties at low temperature T should be the same as for models with more familiar distributions 
(e.g., ± J or Gaussian), but they have the rare advantage that their ground state structure is known. 
The basic assumption of this paper is that a nearest-neighbor Ising spin glass model with any "rea- 
sonable" coupling distribution — i.e., symmetric about zero, and with small or zero weight on 
very large coupling magnitudes — will display qualitatively equivalent thermodynamics at fixed 
d (e.g., presence or absence of a phase transition, number and organization of pure states in the 
spin glass phase, etc.) — at least for very low temperature. While this remains an assumption, it 
is a common one in theoretical spin glass studies — so, for example, the thermodynamics of two 
extreme cases, the ± J and Gaussian distributions, are usually assumed to be the same [TJ7|]. 

All of the models discussed here have the desired properties at finite A (i.e., when T > 0). At 
infinite A (T = 0) they do not, and we make no claims as to whether any ground state properties 
of ordinary spin glasses can be inferred from these highly disordered models. But the interesting 
aspect of the analysis is that, while the thermodynamic properties at nonzero temperature cannot 
be directly solved for, the ground state properties of our models can, and the analysis in Sec. fV| 
enables us to infer properties of realistic spin glass models at low but nonzero temperatures. 

In general, the number of pure states for a given system remains the same or increases (at 
a phase transition point) as temperature decreases. Of course, this general tendency is violated 
at a first-order phase transition, due to phase coexistence. One well-known example of such a 
violation is the g-state Potts ferromagnet with q sufficiently large (depending on d), which at T c 
has its paramagnetic pure phase coexisting with its q ordered pure phases, while below T c the 
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paramagnetic phase is unstable (for rigorous proofs, see [|38|-p0|]). However, it remains the case, 
for these and other systems with first-order transitions, that the number of pure phases above the 
transition is no larger than the number below. More interesting re-entrant behavior occurs in other 
systems — see, e.g., Ref. (A different type of re-entrant behavior occurs in some spin glasses 
in the temperature-concentration phase diagram but this doesn't appear to violate the general 
rule of number of pure states not decreasing as temperature decreases.) 

However, we are unaware of any natural examples where the number of ground states is smaller 
than the number of pure states at arbitrarily low temperature [|42|]. For the spin glass, we are 
looking at an even weaker claim — namely, that the number of pure states in the metastate does 
not decrease as temperature is lowered. That is, it is sufficient to consider only incongruent pure 
states, as discussed in Sec. [53 , 

These considerations enable us to draw a set of interesting conclusions from our analysis of 
these models. The most natural, and obvious, one is that the number of ground states in the zero- 
temperature highly disordered metastate (e.g., with free or periodic boundary conditions) provides 
an upper bound to the number of pure states at very low temperature seen in realistic spin glass 
models. This ground state structure is known: the T = metastate of the highly disordered model 
is supported on a single pair of ground states below eight dimensions, and (uncountably) infinitely 
many above eight |JT^JT7|] . Our analysis therefore provides evidence in favor of the existence of 
no more than a single pair of pure states at low temperatures in realistic spin glass models below 
eight dimensions [|43|]. (This conclusion is consistent, of course, with there being only a single 
pure state, either paramagnetic or otherwise, at all nonzero temperatures in some or all of these 
dimensions.) Our analysis does not allow us to draw conclusions about what happens above eight 
dimensions. 

There are two other logical possibilities, either of which would also be quite interesting. It 
could be that our general intuition about the behavior of the metastate as a function of temperature 
is violated here, so that the number of ground states is smaller than the number of pure states 
in some dimensions. In other words, in these models there might be a jump in the number of 
thermodynamic states, from larger to smaller, in the metastate at zero temperature. We cannot rule 
out this possibility, other to note that the discovery of a class of models in which this occurs raises 
the interesting possibility that it might occur elsewhere also. 

Indeed, a jump presumably does occur in the ordinary 2D EA Ising spin glass — but it goes in 
the other, more natural, direction, i.e. from smaller to larger as temperature is lowered to zero. The 
2D spin glass is believed to be paramagnetic at all nonzero temperatures, and has at least a pair of 
ground states at zero temperature. Our conclusion is that, while we cannot logically exclude the 
possibility of a lowering of the number of thermodynamic states as temperature goes to zero, it 
appears to be less likely than the number remaining the same or increasing. 

We mention in passing a third possibility — that the number of pure states in a typical spin 
glass metastate does not behave in a continuous or even monotonic fashion at all as temperature 
is lowered pT|]. This possibility presents a picture of the low-temperature spin glass phase far 
different from any that have appeared so far in the literature. An extreme version of this possibility 
(bearing some similarity to eigenvalue dependence on disorder realizations in low dimensional 
localization), which we present primarily for illustrative purposes, is as follows. At any (low) 
temperature not depending on the coupling realization J , there would be no broken symmetry 
(and a unique infinite volume Gibbs state) for almost every J, but nevertheless, if a typical J 
were picked first and then the temperature T were varied, there would be a (countably infinite) 
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dense set of temperatures, depending on J, with broken symmetry pure phases for that J at those 
temperatures! However, because there is no evidence of such a picture to date, we do not pursue it 
further here. 
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